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1. Introduction 

In the last two years the several higher-loop calculations of the different finite-length oper- 
ators in 7/ = 4 supersymmetric Yang-Mills (SYM) theory were performed P1]-[|T0|]. These 
calculations usually consist of two part. The first part comes from the Bethe Ansatz for 
the integrable spin chain, which was found first [0]^ in the leading order for the sim- 



plest BMN-operators fllTQ. Generalization to the higher orders together with the investi- 
gations of the integrable structures from the superstring theory side, started in Ref. [^], 
allowed to formulate all-loop Asymptotic Bethe Ansatz (ABA) []l9|]-[36]. However, as was 
shown explicitly in Ref.p7|] by comparison with the predictions from the Balitsky-Fadin- 
Kuraev-Lipatov (BFKL) equation the ABA gives incomplete result starting 

from the four loops for the twist-2 operators and the wrapping corrections should be take 
into account. Thus, the second part is the calculation of the wrapping corrections for the 
finite-length operators. The first such calculation was performed perturbatively for the 
four-loop Konishi operator ^ 0] by taking only relevant Feynman diagrams, following 
Ref. [PI . Then, the wrapping correction for the four-loop Konishi was calculated using 



the generalization of the Liischer formulae |43|] for the finite-size corrections [g]. 

The result of both methods was confirmed with the full direct perturbative calcula- 
tions without any assumptions at all (even without ABA) by evaluation more then 130.000 
diagrams This result is an excellent perturbative test both for the ABA with the dress- 
ing factor and for the different methods of calculations of the wrapping corrections. 

Then, the general form of the four-loop anomalous dimension for the twist-2 oper- 
ators was obtained and the general structure for the leading finite-size correction was 



^Earlier, the similar integrability was opened in quantum chromodynamics in the Regge limit [ p^ |lj 
and for some of operators [p3[p^. 
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understudied [^]^. This result was used to calculate the general form of the five-loop 
anomalous dimension of twist-3 operators which was confirmed by perturbative cal- 
culation in the most simple case 1^. 

At the same time, few models for the complete spectral equations were proposed in 
the form of the Y-system [ i4 , ^] and the Thermodynamic Bethe Ansatz (TBA) equations 
for the ground state p6|]-[p9[]. They give the same results, which were obtained earlier 
with Liischer's formalism. Nevertheless, it is necessary to check these all-loop proposals 
at the higher orders of perturbation theory. New test for the proposed equations was ob- 
tained after the calculation of the five-loop Konishi with Liischer formalism In our 
previous paper with Tomasz Lukowski and Adam Rej [ |lO| ] we have generalized this re- 
sult for the twist-2 operators with the arbitrary number of the covariant derivatives. The 
importance of the obtained result is related with the exceptional properties of the twist- 
2 operators, namely, there is an excellent test for this result, originating from the BFKL 
equation, which gives prediction at any loop order for the analytical continuation of the 
complete anomalous dimension. We have obtained full agreement with both the lead- 
ing order and the next-to-leading order [5C, ^] BFKL predictions. The most important 
consequence of the correctness of obtained result is the general structure of the modifica- 
tion of ABA due to the wrapping corrections. Recently, using the same conditions for the 
ABA modification the result for the five-loop Konishi was computed from the TBA equa- 
tions [52, 53]. Having in hands all necessary formulae and well working programs for the 
calculations of the first Liischer correction at the next-to-leading order for the twist-2 op- 
erators it is naturally to extend them to the analogies calculations for the twist-3 operators, 
as was done earlier at the leading order in Ref.[^]. 

So, in this paper we generalize our previous calculations to the twist-3 operators at 
six-loop order. First, we calculate the part of anomalous dimension, which comes from 
ABA. Then, we write down explicitly all equations that need to be modified at transition 
from the twist-2 case to the twist-3 case. After that, we present the result of finite-size 
corrections with some details of calculations. In the end we will check obtained results 
with known constraints. 



2. The six-loop anomalous dimension from Bethe ansatz 

Twist-three operators belong to the sl(2) sub-sector of the full theory. They can be repre- 
sented by an insertion of M covariant derivatives V into the protected half-BPS state 

Ti {V""^ ZV"^ ZV"""" Z) + ... , mi + m2 + m3 = M. (2.1) 

In the spin chain picture such operators are identified with the states of the non-compact 
s[(2) spin = — I length-three Heisenberg magnet with M excitations. The anomalous 
dimension 7(5) 

7(5) = E72.5^^ 5^ = ^. (2-2) 
^The leading transcendental contribution was obtained earlier by direct perturbative calculation in 
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may be determined to the four-loop order 0{g^) with the help of the ABA [||, and 
can be found in Refs.|37, 54]. Starting from the five-loop order ABA gives the incomplete 
result due to wrapping corrections, which we will discuss in the next section. So, the 
anomalous dimension can be split into the ABA part and the contribution of the wrapping 
interactions 

7(M) = 7^^^(M) + A^(M). (2.3) 

The long-range asymptotic Bethe equations for the s[(2) operators can be found di- 
rectly from the full set of the asymptotic Bethe equations proposed in 1 28, 36] 
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There are M equations for A; = 1, . . . , M which need to be solved for the Bethe roots Uk. 
The variables x^ are related to Uk through 



= u± 
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The function 6{u, v) is the dressing phase and has been conjectured in ||3^]. To the sixth 
order in perturbation theory it has the following form 

e{uk,Uj) = (4C(3)/ -40C(5)/ +420C(7)5'°) iq2{uk) q3{u,) - qsiuk) q2{u,)) 



(2.6) 



- 8C(5)5 iq2{uk) qsiuj) - qbiuk) q2{uj)) 

+ 2Aa5)g^°{qz{uk)qA{uj)-qAiuk)q3iu,))+Oig^^) 



where qr{u) are the eigenvalues of the conserved magnon charges, see |28]. The asymp- 
totic all-loop anomalous dimension is given by 



M 



'i9) = 2g^j: ■ 



(2.7) 



fc=i 



At one loop the Bethe roots Uk are given by zeros of the Wilson polynomial [ p^ , 54 ] 



Pm{u) = 4F3 
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(2.8) 
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1, 1, 1 

Closed expressions for the corrections to the Bethe roots to the three-loop order have been 



obtained in pq ] from the Baxter approach [p6[], which were also extended to the higher- 
loop orders [p7[]. 

In order to obtain closed expressions for the anomalous dimension we solve Eq.(^!^) 
perturbatively for fixed values of the spin M and match the coefficients in an appropriate 
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ansatz, which assumes the maximal transcendentality principle [|^]"^. The basis for the 
ansatz is formed from the harmonic sums, which are defined by the following recurrent 
procedure (see [|63|]) 

S.lM) = f:^S^, ,2.9) 

s.. ..(M) = E^5^^« -(;)■ P-io) 

J=l J 

In the twist-three case the harmonic sums entering into basis should have all positive 
indices and argument M/2 (see Refs.[^ ||]). The number of such sums at the /s-loop 
order for the transcendentality level 2/s — 1 is equal to 4*'"-^. Up to four-loop order we 
need 4^ = 64 harmonic sums and the result for the anomalous dimension can be easily 
foimd (see Refs.[^, 54]). However, the number of sums grows rapidly and already at five 
loops we have 256 harmonic sums in the basis, while at six loops the basis will contain 
1024 harmonic sums. Fortunately, the reciprocity ^A, 65] enters into the game. The pres- 
ence of some structures in the anomalous dimension of twist-2 operators in .A/" = 4 SYM 
theory can be seen at the three-loop order (see Ref.[501). The origin of such structures is 
related with the generalized Gribov-Lipatov reciprocity ^4\ , ^ ]. The reciprocity function 
P^^^(JV) [HH, H] is defined as 

7^^^(M) = P^"'' (^M + ^7^«^(M) j (2.11) 

and considerably simplifies the structure of the anomalous dimension due to additional 



symmetry. Upon substituting the perturbative expansion ( |Z2D , one finds 



P^'^^iM) = f2 5^' 'P2i{M) . (2.12) 



1=1 



The reciprocity-respecting basis appears from the constraints on the harmonic polyloga- 
rithms entering into the evaluation kernel and their Mellin transformations give the nec- 
essary combinations of the harmonic sums []65|,|67|]. These reciprocity-respecting sums are 
defined recursively in the following way (see Ref.[68]): 



Qa = Sa, na,b = i^ai^b) , ^a,b,c = ^ai^b.c) , (2-13) 

Wa(5'i,c) = Sa,b,c - ^SaAb,c , u A m = sign(n) sign(m) (|n| + |m|) . (2.14) 



Several theorems were proved for the corresponding combinations of harmonic sums [ |65t 



67, ^ which allowed to reduce considerably the number of harmonic sums entering 



^The hypothesis about maximal transcendentality principle [ pSP was confirmed by direct perturbative 
calculations at the two-loop order [^] and then successfully applied at the three-loop order [|60|], when corre- 



sponding results were obtained in QCD ^ ||]. 
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into the basis. In this way, the result for the five-loop anomalous dimension for the twist- 
three operators was obtained in Ref.[01. 

The same basis can be constructed in a more elegant and natural way with the help 
of binomial harmonic sums. These sums come both from the perturbative calculations 
of the anomalous dimensions of twist-2 operators [61, and from the solution of the 
Baxter equation for the corresponding spin chain [|69|,|5^]^. In our previous paper [ [T0| ] we 
found, that the basis from the binomial harmonic sums is equivalent to the reciprocity- 



respecting basis [ 6£ ]^ and we used the binomial harmonic sums for the calculation of five- 
loop anomalous dimension of the twist-2 operators both for ABA and for the wrapping 
correction parts. We define the binomial harmonic sums through (see |]63[]) 



(i\r) = (-i)^^(-i) 




(2.15) 



where 5ii_..._i^ are the nested harmonic sums defined in ( 2.1C| ). One of the interesting fea- 
ture of these sums is that they are defined only for positive values of the indices ii,. . . ,ik- 
More interesting, that there are binomial harmonic sums, which are expressed through 
the usual harmonic sums with all positive indices, for example 



52,1,2,2,1,2,1 



3,5,3 



85; 



3,5,3 



4 5: 



3,8 



4 58,3 + 25 



11 



(2.16) 



Namely these sums will form the basis for reciprocity functions Vz of ABA part for twist-3 
anomalous dimension! The number of such sums, which will form the basis, for tran- 
scendentality level k = 1, 3, 5, 7, 9, 11 are equal to 1, 2, 5, 13, 34 and 89 respectively^, 
instead of 4*'"-^ as for the usual harmonic sums with all positive indices. Thus, the basis is 
significantly reduced and the reciprocity functions Vi can be easily calculated. 

From Eqs.(^lT|) and (|1|) one can find, that at six-loop order the reciprocity function 
Pi2 is related with the anomalous dimension 712 (see Appendix B of Ref.[|^ for five loops): 



7^12 (M) = V12 



-T-jrational 
'12 



(2.17) 



-T-jrational 
'12 
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'^12 
^12 



rational 
7l2 



76^ + 2747r*'°""' + 2 727i^g' 



..rationaA ' 



+ - 



7I + 6 727674 + 37|7r*"°"')" 



96 

--^71(371 +2 7276 



768 
7,?' 

71?' 
7,?' 



+ 



6144 



7I74 



737280 



7.7s«='+7.7jf)' + ^ 
7.72")' . 



7l7l<")" 



(2.18) 

(2.19) 

(2.20) 
(2.21) 



We will give more details about this in the forthcoming paper. 

^Relations between the binomial and the nested harmonic sums together with relations 
between the binomial and the reciprocity-respecting harmonic sums can be found under 

http : //thd . pnpi . spb . ru/ ~ velizh/ Sloop / and http : / /thd . pnpi . spb .ru/~velizh/51oop3/ 



These numbers are the subset of the Fibonacci numbers: 0, 1, 1, 2, 3, 5, 8, 13, 21, 34, 55, 89... 
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where each prime marks derivative over M. 

We have found the following result for Vi2{M) 

-2532 §2,2,2,2,2,1 + 2160Si^i_2,2,2,2,l + 1408 §1^2,1,2,2,2,1 + 1056 §1^2,2,1,2,2,1 



Trjrational 
^^12 



+ 1264 §1 

+944 §2,1 
+ 1136 §2 
-1376§i 
-736 §1,1 
-960 §1,1 
-960 §1,2 
-320 §1,2 
-256 §2,1 
-720 §2,1 
-864 §2,1 
-864 §2,2 
+384 §1,1 
+384 §1,1 
+ 192 §1,1 
+256 §1,2 
+448 §2,1 
-512 §1,1 



2,2,2,1,2,1 + 1856 §1,2,2,2,2,1,1 + 1248 §2,1,1,2,2,2,1 + 768 §2,1,2,1,2,2,1 
,2,2,1,2,1 + 1280 §2,1,2,2,2,1,1 + 784 §2,2,1,1,2,2,1 + 832 §2,2,1,2,1,2,1 
2,1,2,2,1,1 + 864 §2,2,2,1,1,2,1 + 1 136 §2,2,2,1,2,1,1 + 1504 §2,2,2,2,1,1,1 
1,1,1,2,2,2,1 — 896 §1,1,1,2,1,2,2,1 — 1120 §1,1,1,2,2,1,2,1 — 1312 §1,1,1,2,2,2,1,1 
2,1,1,2,2,1 — 960 §1,1,2,1,2,1,2,1 — 1088 §1,1,2,1,2,2,1,1 — 640 §1,1,2,2,1,1,2,1 
2,2,1,2,1,1 — 1088 §1,1,2,2,2,1,1,1 — 512 §1,2,1,1,1,2,2,1 — 864 §1,2,1,1,2,1,2,1 
1,1,2,2,1,1 - 656 §1,2,1,2,1,1,2,1 - 848 §1,2,1,2,1,2,1,1 " 1056 §1,2,1,2,2,1,1,1 

768 §1,2,2,1,1,2,1,1 - 992 §1 2.2.1.2.1.1.1 - 768 §1 



1,1,1,2,2,1 



800 § 



2,1,1,1,2,1,2,1 



704 § 



2,1,1,1,2,2,1,1 



624 § 



2,1,1,2,1,1,2,1 



1,2,1,2,1,1 - 928 §2,1,1,2,2,1,1,1 - 512 §2,1,2,1,1,1,2,1 - 672 §2,1,2,1,1,2,1,1 
2,1,2,1,1,1 - 864§2,l,2,2,l, 1,1,1 - 96 §2,2,1,1,1,1,2,1 " 512 §2,2,1,1,1,2,1,1 
1,1,2,1,1,1 - 864 §2,2,1,2,1,1,1,1 - 640 §2,2,2,1,1,1,1,1 + 1024 §i,i,i,i,i,2,l,2,l 
1,1,1,2,2,1,1 + 640 §1,1,1,1,2,1,1,2,1 + 768 §1,1,1,1,2,1,2,1,1 + 512 §1,1,1,1,2,2,1,1,1 



1,2,1,1,1,2,1 + 640 §1,1,1,2,1,1,2,1,1 + 832 §1,1,1,2,1,2,1,1 

2,1,1,1,1,2,1 + 512 §1,1,2,1,1,1,2,1,1 + 832 §1,1,2,1,1,2,1,1 

1,1,1,1,2,1,1 + 704 §1,2,1,1,1,2,1,1,1 + 832 §1,2,1,1,2,1,1,1 

1,1,1,2,1,1,1 + 704 §2,1,1,1,2,1,1,1,1 + 640 §2,1,1,2,1,1,1,1 
1,1,1,2,1,1,1,1 • 



1 + 256 §1,1,1,2,2,1,1,1,1 
1 + 832 §1,1,2,1,2,1,1,1,1 
1 + 640 §1,2,1,2,1,1,1,1,1 
1 + 640 §2,1,2,1,1,1,1,1,1 
(2.22) 



The functions 'Pi2^\ 1^12^"^ and Vi2''^ are easy to determine 

7-'l2^^= — 32(45§i,2,2,2,l+ 21 §2,1,2,2,1+ 19§2,2, 1,2,1+ 39 §2,2,2,1,1 — 6 §2,1,1,1,2,1 — 8 §2,1,1,2,1,1 
— 32 §1,1,2,2,1,1 — 10 §1,2,1,1,2,1 — 12 §1,2,1,2,1,1 — 28 §1,2,2,1,1,1 — 24 §1,1,1,2,2, 1 

-20 §1,1,2,1,2,1- 16 §2,1,2,1,1,1- 12§2,2,i,i,i,i+24§i,i,i,2,i,i,i+24§i,i,2,i,i,i,i), (2.23) 

T'fr = -64 (29 §1,2,2,1 + 5 §2,1,2,1 + 28 §2,2,1,1 - 20 §1,1,2,1,1 - 19 §1,2,1,1,1 + 3 §2,1,1,1,1) , (2.24) 

= -3360 (§1,2,1 + §2,1,1) . (2.25) 

The expression for the six-loop anomalous dimension of twist-three operators from ABA 
in the canonical basis of the nested harmonic sums ( p.lO ) can be found on our web-page. 

3. Wrapping correction 

In this section we will write explicitly all equations, which should be modified in Ref.[lC] 
for the calculations in the case of twist-three operators. The wrapping correction is cal- 
culated by evaluating the first Liischer correction at weak-coupling along the lines ad- 
vocated in ^ 1^. The interested reader should read this subsection together with our 
previous paper [[To[], where all definitions and a lot of details can be found. 
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First of all, the length of twist-three operator is equal to L = 3, which according to 
general formulae for F-term (see Subsection 5.2 in Ref. ^ for details) 

L 



E 



dq ( z~ 
,27r \z- 



^(-1)^' [5g_i(g,u,)5g_i(g,u,,)]J[;;| (3.1) 



renders the exponential part to be of the form 



645^ 



(g2 + Q2)3 



24 



This leads to the following modification of the last factor in Eq.(5.6) of [ [10[ | 

16 64 



— j> 



and 



(g2 + Q2)2 (g2 + g2)3 



C'°^(M,?) = 645i2^^^^'^^ 



64 



RM{q,Q) (92 + Q2)3 



(3.2) 



(3.3) 



(3.4) 



where Yq'^''^\m, q) is the integrand for the five-loop wrapping corrections to the anoma- 
lous dimension of twist-3 operators [0] 



aL^°)(m) 



E 

0=1 



(3.5) 



with Tm(?, Q) and i?M(9, Q) are given by 



i?ff(g,g) = P 



(0) /g-i(-l + Q)\ ^(0) fq + i{-l + Q)\ <o) fq-i{l + Q) 



M 



M 



M 



XP; 



(0) fq + i{l + Q) 



M 



(-1) 



M 



,(0)/ Q-i{Q- 1) 



2j-iq-Q 2{j + l)-iq-Q^ 



+ ij 



(3.6) 
(3.7) 



Farther, the one-loop energy of twist-three operators differs from the twist-two one. It is 
given by 

85i(f) 



M 



(3.8) 



and we should replace everywhere in 5i(M) argument M to M/2. Moreover, the one-loop 
Baxter fimction Eq.(5.9) of [|T0|], for twist-three case is given by 



_M M,^ 1,- 1 
2 1 P -^1 2 '"■1 2 



2 

1, 1, 1 



lU 



(3.9) 



The integrand y(^2,o) £qj. ^.j-^g six-loop wrapping corrections to the anomalous dimen- 
sion of twist-3 operators can be split by analogy with the five loops for the twist-2 oper- 
ators into the matrix part, scalar part, the exponential term and the dressing factor (see 
details in Refs.g and 
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For the scalar part of y(i2,o) is easy to find, that 5_2(M) term in Eq.(5.11) of 0] 
drop out (and do not forget to change the argument of Si in the last term) 



5gL(M,g,Q) 

,(0), 



4? /pW'(^i4M) p(0)'(2^,(i±Ql)\ I6g5i(f) 



The exponential part Eq.(5.15) of [ p^ ] should be modified according to Eq.(^) as 

T(6)(g,Q) 24 



(3.10) 



(3.11) 



The two-loop solution P)^'{u) Eq.(5.25) of [jlOj] can be taken for twist-3 from Ref.[E 



25i(f ) (45i(f ) - 25i(M)) + 3 52(f) ] P^°\u) 



■f , f + i + 2<5, l+m, I 
1 + (5, 1 + 6, 1 



5=0 



3 

2 ds^ 



1P3 



2 ' 2 ^ ' 2 ^ 2 



(3.12) 



5=0 



1 + ^, 1-5, 1 

In the modification of ABA part it is necessary to change the powers in Eq.(5.32) of [ |10| ] 

3 M 



BY{uk) 



n 



Uk 



Uj + I 



(3.13) 



and in Eq.(5.33) of ^ 



= £ ^ (it) str|5Q(g, ui) . . . a5Q(g, Uk) . . . Sqiq, um)^ . (3.14) 

We have checked, that the six-loop integrand obtained with formulae from Section 
5 of Ref.[10], modified according to this section, gives zero for the sum over dynamical 
poles in M = 2 and M = 4 cases''. 

3.1 Calculation and result 

For calculations we use the same method as in Ref.[|^]. With formulae from Section 5 of 
Ref.[p^, modified according to the previous section, we obtain the integrand for given 
M, put Q = 1, perform numerical integration over q and repeat this procedure for Q = 
2,3,.... Then, the high-precision result is used to find the coefficients in the following 
transcendental basis 



{as), c(3)^ c(5) as), a7), c(3)^ as), as), 



(3.15) 



using EZ-Face [|7^] and/or PSLQ algorithm [|7ip implemented on MATHEMATICA [|7|] 
(EZ-Face is limited with 100 digits of accuracy). This allowed us to conjecture the form 



thank Tomasz Lukowski and Adam Rej for the check of the cancelation. 
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of the functions multiplying zeta functions. For higher values of M the non-rational part 
can be subtracted from the full result and it is sufficient to perform computations to the 
accuracy that allows to rationalize the result. In this way we obtain 25 values up to M = 
50. 

For the reconstruction of the full M-dependence we assume that the wrapping cor- 
rections preserve the reciprocity symmetry. This implies that a part of the wrapping cor- 
rection may be found from the lower order results 

Al") = ^(AL^%)' + Pr2. (3.16) 



2 

where 72 is the one-loop anomalous dimension and A^°^ is the five-loop wrapping cor- 
rection found in [^] (all harmonic sums being of argument M/2) 

72 (M) = V2iM) = 8 5i(M/2) =4Si(M/2), (3.17) 

aL^°^ = P^o = -645^° Sf{35Ci7) - 4052C(5) + (-854 + 1652,2)C(3) 

+257 - 452,5 - 253,4 - 454,3 - 256,1 + 852,2,3 + 453,3,1) (3.18) 

= -64 (35 C(7) + 4 f^3,3,i + 8 ^^2,2,3 + 24 C(3) Q2,2 - ^7) • (3.19) 



Here Qm are the complementary reciprocity-respecting harmonic sums defined through 
complementary harmonic sums 5m as (see Appendix A in Ref.[|^] and Ref.[||]) 

5a,m(M) = 5a,m(M) - 5a(M)5m( 00) . (3.20) 

We suggest, that the remaining part V12 should have the same general structure as for 
twist-two case [ M 



7^1-2 = 2V^T + 2V2[2Vi + —7^2 J (35C(7) - 4052^(5) + (-854 + 1652,2)C(3) 



+257 - 452,5 - 253,4 - 454,3 - 256,1 + 852,2,3 + 453,3,iJ , (3.21) 
with T and are given by 

r = c(9) rc(9) + c(3)' r^(3)3 + C(5) c(3) rc(5)c(3) + C(7) rc(7) + a^f 

+C(5) rc(5) + C(3) rc(3) + T^ational , (3.22) 

Vi = -853 . (3.23) 

Applying the principle of maximal transcendentality [ |5^ ] we conclude that the tran- 
scendentality of the components 7^(9), 7^(3)3, 7^(5)C(3)/ 7^(7)/ 7^(3)=/ 7^(5), 7^(3)/ 7;ationai 
should be equal to 2, 2, 3, 4, 5, 6, 8 and 11 respectively. The lowest-transcendentality 
functions, 7^(9), 7^(3)3/ 7^(5)^(3)/ 7^(7) and 7^(3)2 may be obtained from first few values 

(3.24) 
(3.25) 
(3.26) 
(3.27) 
(3.28) 





= 64512 Sl , 


Tc(3)3 


= 0, 


'^C(5)C(3) 


= -7680 Sf , 




= -62720 5^52 




= 6144 5^52 . 



-9- 



Unfortunately, as can be seen from Eq.( |3.18 ) the basis for other Ti will contain the 
harmonic sums with even numbers among indices, in contrast to the ABA part. Thus, the 
basis will extended considerably, even if one take into account the reciprocity-respecting 
parity even sums [^. Moreover, the sums are appeared, which depend on M, not on M/2. 
This can be seen by applying MATHEMATICA function Factor Integer to the denomina- 
tors of the obtained values. For example, for M = 8 we have found: 

A<")(8) = 280000fl9) - C,3)f(5) - f(7) + C(3)^ 

9539440625 , 294429940625 , 277950256725625 

+ C(5) C(3) + (3.29) 

40824 ' 2939328 ' 6772211712 ^ ' 

and Factor Integer gives for the denominator of C(5)'s coefficient: 

40824 = 2^ X 3^ X 7^ (3.30) 

For argument of harmonic sums equal to M/2 the last number 7^ should be absent if M = 
8. Thus, the basis expands extremely. Fortunately, there is a nice observation in Appendix 
B of Ref.[^ about relation of the exact result with the functions Pm and Tm(9, Q) from the 
integrand for twist-two case at four loops. We assume that similar properties will remain 
for the operators of twist-3. Following Ref.[^, we rewrite Tm(9, Q) as 



TM{q,Q) = E 



Q-l r 1 1 T a 



2j -iq-Q 2{j + l)-iq-Q 

(|(? - ^(g - 1))) iPM(i(g + i(g-i)) 



PM(-iq-tiQ-l)) + ij 



q-iQ q + iQ 



+ fMiq,Q), (3.31) 



- ^ ^ '5=,iPM(|(?-^(g-i)) + ij)-PM(i(g-i(g + i)) + u, ^ ^ 

expand Tniq, Q) over Q for a lot of M and reconstruct the general dependence on M as a 
combination of harmonic sums for each power of Q. As can be seen from Appendix^ for 
example, in the case of Q-'- the harmonic sums 5_2,i(M) and 52,i(M) enter into answer. 
Thus, to reconstruct the full M-dependence for C(5)'s term we write down the following 
basis (if not written explicitly, the argument of harmonic sums is M/2): 

Ql{Qi, Q3,u^2,2, ^1,3, ^^2,i,i> f^i,2,i, f^if^-2,i(M), QiQ2,i{M)} (3.33) 

and we have found 

rc(5) = -512 Ql (S Qi (8 ^^2,1 + ^s) - 50 ^^2,2 - 6 ^^3,1 + 20 ^2,1,1 + 15 ^^4) 

+40960 QI (n_2,i(M) + n2,i(^)) • (3-34) 



For 7^(3) part we add to the standard basis from the harmonic sums with argument M/2 

Ql {Qq, ^5,1, ^4,2, f^3,3, f^2,4, ^1,5, ^^4,1,1,^^1,4,1,^^3,2,1, ^^3,1,2, f^2,3,l, ^^2,1,3, ^^1,3,2, f^l,2,3, 
^^2,2,2, ^^3,1,1,1, ^^1,3, 1,1, f^2,2, 1,1, f^2, 1,2,1, f^2, 1,1,2, ^^1,2,2,1 , ^^1,2, 1,2} (3.35) 
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the following combinations of harmonic sums with argument M 



^2,2,l{M) = fi2,2,l(M) + Q-2,2AM) + ^2,-2,l{M) + f^-2,-2,l(M) , (3.36) 
f22,l,2(M) = n2,i,2(M) + n_2,i,2(M) + n2,i,-2(M) + n_2,l,-2(M) (3.37) 

multiplied by fif (M/2). With this basis we have obtained 

Tc(3) = 1024^^(^1(8(^2,1,2+ 3 f22,2,l) -^^5) - f^5,l + 2^2,4 + 3^4,2 +4^2,1,1,2 

+4 ^2,1,2,1 + 8 02,2,1,0 - 32768 fi? (n2,i,2(M) + 3 fi2,2,i(M)) . (3.38) 



To reconstruct the rational part we need the basis with at least 83 reciprocity-respecting 
harmonic sums with argument M/2 and at least two following combinations of the reci- 
procity-respecting harmonic sums with argument M (see details in Appendix^: 

^^2,2,l,3(M) = fi2,2,l,3(M) + n_2,2,l,3(M) + n2,-2,l,3(M) + n2,2,l,-3(M) + n_2,-2,l,3 (M) 
+^^2,-2,l,-3(M) + n-2,2,l,-3(M) + n_2,-2,l,-3(M) , (3.39) 

^^2,1,2,3(M) = n2,l,2,3(M) + n_2,l,2,3(M) + fi2,l,-2,3(M) + n2,l,2,-3(M) + n_2,l,-2,3 (M) 
+f^2,l,-2,-3(M) + f^-2,l,2,-3(M) + f^-2,l,-2,-3(M) . (3.40) 

So, we should calculate 85 values up to M = 170. This is far beyond computational capa- 
bilities. However, as we well know from the previous results, coefficients in the front of 
the reciprocity-respecting harmonic sums should be integers with a lot of zeros among 
them. We can try to find these coefficients with the Number Theory. The most sim- 
ple algorithm for the solution of given problem is the LLL-algorithm []73|], which is re- 
alized inMATHEMATICA with function LatticeReduce. We calculate the values of all 85 



reciprocity-respecting and special sums in the basis [A.7[]A.9| up to M = 48. So, we have 
24 equations in the linear system for 85 variables together with our corresponding results 
with subtracted low order anomalous dimension according to Eq.( p.l6| ). We eliminate 23 
variables and we remain with one equation on 62 variables. According to the realization 
of LLL-algorithm^ we replace the last column in 64 x 64 unity matrix with the last equa- 
tion of our system (62 + 1 numbers) and zero as the last element. After four minutes^ of 
calculations MATHEMATICA gives: 

T^ational = —512 ^i (^i (^2,6 + ^^5,3 — 8(f^2, 1,2,3 + ^2,2,1,3)) + + f^2,l,6 + ^2,6,1 — ^3,3,3 
-f^3,5,l + ^5,1,3 - f^5,3,l + ^7,1,1 " 2 f^2,2,5 " 2 f^2,4,3 " 3 ^4,2,3 " 4 f^2, 1,1,2,3 
— 4^^2,1,2,1,3 — 4 f22, 1,2,3,1 — 4^2,2,1,3,1 — 4^23,1,1,3,1 — 4f^3,3,l,l,l — 8 f^2,2,l,l,3 



-8 f^2,2,3,i,i j - 65536 ni (C!2,i,2,3(M) + fi2,2,i,3(M) j . (3.41) 

Note, that LatticeReduce gives 64 possible solutions but all other solutions contain 
large integers and only few zeros. We have checked the obtained result for M = 25 and 



^see Application on |http ://reference. wolfram . com/ mat hematica/ref /LatticeReduce . html" 



The computation of M = 50 requires about 200 hours. 
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have found the full agreement, which can serve as confirmation of the correctness for 
obtained result. 

Now we are ready go to the complementary functions Qa for testing of T on the parity 
following Ref.[0|. We have found 



C(3)3 



64512 Ql , 
0, 



'^C(5)C(3) 



-7680 Q 



1 1 



rc(5) = -1536 (5 f^if^s - 19 - ^3,1) , 



(3.42) 
(3.43) 
(3.44) 
(3.45) 
(3.46) 

(3.47) 



where we don't write the expressions for 7^(3) and T^ationai as the harmonic sums, entering 



into Eqs.( |3.38| ) and ( |3.4lD have only odd numbers through indices or two even numbers 
with other odd numbers through indices. According to theorem from Ref.[^] all comple- 
mentary functions in the expressions for are parity-even, then T is parity-even also and 
satisfy reciprocity. With help of the SUMMER [||] and HARMPOL packages for FORM [ |7^ ] 
and the HPL package [|7|] for MATHEMATICA, we checked, that the combinations of har- 



monic sums in the special sums Eqs.( |3.38 ) and ( 3.41 ) translated into polylogarithms satisfy 
Gribov-Lipatov relation [ |65| ] up to ^'s terms. However, the complete analysis of the full 
result for the six-loop anomalous dimension of the twist-3 operators is unavailable with 
these programs due to their limitations. Nevertheless it seems, that the reciprocity will 
work for the full result also. 



4. Large M asymptotic and analytical continuation 



In this section we check our result for the six-loop anomalous dimension of the twist-three 
operators against the known constraints. The M — > 00 limit can be easily calculated with 
use of the SUMMER package [||] for FORM 



M 



lim aL^^^(M) = 0. 



(4.1) 



This means that, the wrapping effects do not influence the scaling function. 

Moreover, there are predictions from the empirically established resummation formu- 
lae for the analytical continuation to M = — 2 + w from Ref.[^ 



7 = -«:^ 




which give at six loops: 



7l2 



8 144^, , 
741+^^(2). 



Q 

t — — 



(4.2) 



(4.3) 



With help of the SUMMER [|63|] and HARMPOL [|7^] packages for FORM [|7|] and the HPL 
package [ |76| ] for MATHEMATICA, we have found the following results for the analytical 
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continuation to M = —2 + w from ABA part, from the wrapping corrections part with- 
out special sums and from the special sums (the last terms) in the wrapping corrections 
correspondingly: 

1288 6800 ,,, 768 5120 ,,, 512 1536 ,,, 



Summing up all three pieces we obtain full agreement with (4.3), which can serve as ad- 
ditional test for the correctness of obtained result. 

In the end we want to write explicitly few first values of the six-loop anomalous di- 
mension of the twist-3 operators from our general result ( |2.18D -( ^!25| ), ( |3.16| )-( ^!28|) , ( |3.34| ), 



( p.38D and ( |3.41| ). The anomalous dimension of the most simple operator with M = 2, 



which is the analog of Konishi operator in twist-2 case, is equal to 

7i2(M = 2) = -53016 - 18176 C(3) - 16128^(5) - 13440^(7) 

+3840 + 6144 C(3)^ - 7680 C(5) C(3) + 1024 C(3) - 3584 C(5) - 53760 C(7) 
+64512 C(9), (4.5) 

where the first line is the contribution from ABA. In the sum it gives 

7i2(M = 2) = -49176 + 6144 C(3)^ - 7680 C(5) C(3) - 17152 C(3) - 19712 C(5) 

-67200 C(7) + 64512^(9). (4.6) 

For M = 4 we have found: 

3 1^68363 

7i2(M = 4) = — 36018 C(3) - 29916^(5) - 22680 C(7) 

256 

127431 

+ + 25920 C(3)^ - 25920 C(5) C(3) - 14364 C(3) + 36216 C(5) 

8 

-176400 C(7) + 145152 C(9) . (4.7) 



And for M = 6, which ABA part (divided by 2 ) can be found in Ref.|57], we have: 



2720281112987 317437583 , 3175148 „ , 773080 

7i2(M = 6) = C(3) C(5) C(7) 

' ^ ' 15116544 6561 ' 81 27 ' 

68920603499 4174016 425920 , „ , 1650710908 „ , 
+ 2361960 + ^ ^<=' ^<^' 3280^ ^<^' 

+^«5)-^C(7) + 2ie832«9). ,4.8, 
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A, Special sums 



In this Appendix we will show the appearance of the special harmonic sums in the basis 
for the reconstruction of the general expression for the wrapping corrections. In Appendix 
B of Ref.[^ it is shown, that the results for some of ^'s contributions to the four-loop 
anomalous dimension of the twist-two operators from the wrapping corrections can be 
obtained exactly with the expansion over Q of the functions Pm and TmIq, Q) entering 
into integrand. In our case we will concentrate only on Tm{q, Q) function, which can be 
written according to in the following form 



Q-i 

TMiq,Q) = E 

i=o 



2j-iq-Q 2{j + l)-iq-Q 



PM(^{q-i{Q-l)) + ij 



iPm Uiq - iiQ - m iPM U(g + «(g-l)) 



q-iQ q + iQ 

Pm {^{q - i{Q - 1)) + ij) - Pm [\{q - i{Q + 1)) + ij 

i=i 



+ fM(g,g), (A.l) 



Tm(., q) = E ^ ^-^-r^ ^ • (A.2) 



We substitute two-loop Baxter function P^ from Eq. (t3.12| ) instead of Pm, expand Tm(5, Q) 



over Q for a lot of M and for each power of Q reconstruct the combination of harmonic 
sums, which give expression for arbitrary M. For Q = the transcendentality level for ob- 
tained result will be the same, as transcendentality level of two-loop Baxter function, that 
is equal to 3 and we can easily found (argument of harmonic sums is M/2 if not written 
explicitly) 

Q° : 653 - 25i (-45i(M)5i + 851^ + S2) . (A.3) 

Similar structures come from the expansion of Pm, so the harmonic sums with argument 
M will drop out in the expression for the anomalous dimension. For the first power of Q 
the transcendentality level for obtained result is 4 and we have found 

: 2 (-125i,3 + 352,2 - 853,i + 75i,i,2 + 35i,2,i - 52,i,i + 45^^ + 454 

-45i(M)5i (51^ + 52) + 165i (5_2,i(M) + 52,i(M)) . (A.4) 

The sums in the last term are the special one, which do not appear in the expansion of Pm 
and can enter into the anomalous dimension. So, we will add such sums multiplied by 
common factor 5^ (M/2) to the basis for the transcendently level 6, which is relevant for 
the (^(5) contribution. For the second power of Q the result of expansion Tm(5, Q) 

Q2 : ^5i(M) (5f - 5i53) - 25i (253,1 + 54) - 2 (553,2 - ^2,3 + 654,i - 1253,i,i + 55) 
-^Sl + 2S2Sl + '^SiSl (A.5) 
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does not contain of the special sums, while the result for the third power of Q 



: -^5i(M)5i (652,2 + Sf - 4:3^3, - 33^) + 25? (852,2 - 253,i + 54) 
+645i (5_2,i,-2(M) + 5_2,i,2(M) + 52,i,-2(M) + 52,i,2(M)) 



+-5i (-552,3 + 453,2 + 354,1 - 652,1,2 + 652,2,1 - 1253,1,1 + 555) 



-1552,4 + 2053,3 + 954,2 + 1255,1 + 652,2,2 - 1252,3,1 - 3653,1,2 - 3653,2,1 




(A.6) 



contains the special sums again. This observation and its generalization to the higher 
special sums allows us to assume the minimal set of the most simple special sums, which 
can enter into the basis for the rational part of T- So, the minimal basis for the rational part 
of T has 83 reciprocity-respecting sums of the transcendentality level 9 with argument 
M/2 

{Qq, Qlfifi, ^^1,3,5, ^^l,4,4i ^^l,5,3i ^^1,6,2, ^^1,7,1 1 ^^2,1,6 1 ^2,2,5, ^^2,3,4i ^^2,4,3, ^2,5,2, ^^2,6,li 
^3,1,5, ^3,2,4, ^^3,3,3, ^3,4,2, ^^3,5,1 , ^4,1,4, ^^4,2,3, ^4,3,2, ^^4,4,1, ^^5,1,3, ^5,2,2, ^^5,3,lj ^6,1,2, 
^6,2,1, f^7,l,li ^1,2,1,1,4, f^l,2,l,2,3i ^1,2,1,3,2, ^^1,2,1,4,1 , ^^1,2,2,1,3, ^1,2,2,2,2, ^^1,2,2,3,1 1 ^1,2,3,1,2, 
^^1,2,3,2,1, ^^1,2,4,1,1, ^^1,3, 1,1,3, ^^1,3, 1,2,2, ^^1,3, 1,3,1, ^^1,3,2, 1,2, ^^1,3,2,2,1, ^^1,3,3,1,1 , ^^1,4,1, 1,2, 
^1,4,1,2,1, ^^1,4,2,1,1, ^^1,5, 1,1,1, ^2,1,1,1,4, ^^2, 1,1,2,3, ^2,1,1,3,2, ^^2, 1,1,4,1, ^^2, 1,2, 1,3, ^2,1,2,2,2, 
^^2, 1,2,3,1, ^^2, 1,3, 1,2, ^^2, 1,3,2,1, ^^2, 1,4,1, 1 , ^^2,2, 1,1,3, ^^2,2, 1,2,2, ^^2,2, 1,3,1, ^^2,2,2, 1,2, ^^2,2,2,2,1, 
^2,2,3,1,1, ^2,3,1,1,2, ^^2,3, 1,2,1, ^^2,3,2,1,1 , ^^2,4,1,1,1 , ^^3, 1,1, 1,3, ^3,1,1,2,2, ^^3, 1, 1,3, 1 , ^3,1,2,1,2, 
^^3, 1,2,2,1, ^^3, 1,3, 1, 1, ^^3,2, 1,1,2, ^^3,2,1,2,1 , ^^3,2,2,1,1 , ^^3,3,1,1,1 , ^^4, 1,1, 1,2, ^^4,1,1,2,1 , ^^4,1,2,1,1 , 
^^4,2,1,1,1, f^5, 1,1, 1,1} (A.7) 



and two special sums with argument M 

^^2,2,l,3(M) = f^2,2,l,3(M) + f^_2,2,l,3(M) + f^2,-2,l,3 (M) + f^2,2,l,-3(M) + f^_2,-2,l,3 (M) 



^^2,1,2,3(M) = f^2,l,2,3(M) + f^_2,l,2,3(M) + f^2,l,-2,3 (M) + f^2,l,2,-3(M) + ^-2,1,-2,3 (M) 
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